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We calculate the non-resonant Raman response, the single particle spectra and the charge-spin 
configuration for the electron doped t — t' Hubbard model using unrestricted Hartree-Fock calcu- 
lations. We discuss the similarities and differences in the response of homogeneous versus inhomo- 
geneous structures. Metallic antiferromagnetism dominates in a large region of the U — n phase 
diagram but at high values of the on-site interaction and for intermediate doping values inhomoge- 
neous configurations are found with lower energy. This result is in contrast with the case of hole 
doped cuprates where inhomogeneities are found already at very low doping. The inhomogeneities 
found are in-phase stripes compatible with inelastic neutron scattering experiments. They give an 
incoherent background in the Raman response. The B2g signal can show a quasiparticle component 
even when no Fermi surface is found in the nodal direction. 
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I. INTRODUCTION 

Electron and hole doped cuprates share a similar 
phase diagram with antiferromagnetic and superconduct- 
ing phases. However there are quantitative differences 
between them. While electron-doped show commensu- 
rate antiferromagnetismi in a large region of the U — n 
phase diagram there is a quick suppression of antiferro- 
magnetism in hole-doped cuprates. There it is well es- 
tablished that inhomogeneities appear when doping the 
Mott insulator. Among them the stripes, where holes 
are aligned in a row, have acquired importance as a pos- 
sible scenario to explain the phase diagram of the hole 
doped cuprates--. The strongest experimental support 
for stripes comes from neutron scattering experiments^ 
that show incommensurate peaks. They have been inter- 
preted as coming from anti-phase stripes where the line 
of holes separate antiferromagnetic domains of opposite 
sign. The case of the electron doped cuprates is less clear. 
The signal obtained in neutron scattering experiments is 
commensurate^ what would be compatible with in-phase 
inhomogeneities. This possibility can convene with theo- 
retical predictions that incommensurate solutions will be 
more favored in hole than in electron dopingSiSiLSiS. Very 
recent measurements of thermal conductivity and neu- 
tron scatteringiS in electron doped cuprates have been 
interpreted as in-phase stripes while experiments of Nu- 
clear Magnetic Resonanceii were seen as signatures of 
inhomogeneous structures. A low temperature pseudo- 
gap detected with tunneling spectroscopjii^ points to the 
presence of a ground state other than superconductivity 
whose origin could rely in inhomogeneous configurations. 

In this work we perform a systematic analysis of homo- 
geneous versus inhomogeneous structures in the electron 
doped t — t' Hubbard model using unrestricted Hartree- 
Fock (UHF) calculations. Stripes and inhomogeneities 
were predicted as a solution of the Hubbard model with 
t' = using UHF approximation^^. Within this frame- 
work we compute the spectral function and the non- 
resonant Raman response of the different configurations. 
In the study of commensurate antiferromagnetism it is 



interesting to have charge sensitive probes as the non- 
resonant Raman response which can detect charge or- 
der due to the possible selection of different regions of 
the Brillouin zonei^. Inhomogeneous solutions in elec- 
tron doped have been studied previously with UHF in 
the three band Hubbard modeli^ where they found that 
diagonal anti-phase stripes were the configuration lowest 
in energy. So far there are no evidence of these stripes in 
the experiments. We will use the simplest t — t' Hubbard 
model which has shown to be able to explain the pseudo- 
gap and the phase diagram of electron dopcd^^ cuprates 
and to give the correct low energy physics when com- 
pared with the two-band Hubbard modeliS. The t — t' 
Hubbard model has also been studied with UHF in the 
past22iSi where it was pointed out that t' favors homo- 
geneous solutions in the electron doped case in contrast 
to the hole-doped. This was claimed to be due to the 
fact that t' does not frustrate the antiferromagnetism in 
the electron doped case as was also proven with other 
techniques as exact diagonalization in the electron doped 
t — t' — J modelMi^ and in Quantum Monte Carlo cal- 
culations of the t — t' Hubbard model2iS4. Metallic an- 
tiferromagnetism has been found in experimental. Our 
calculations agree with the former results that homoge- 
neous antiferromagnetism is more stable in the electron 
than in hole doped case but in this systematic study 
we have found that for large enough values of the on- 
site interaction U = 6t ~ 8t and at intermediate dop- 
ing in-phase inhomogeneous configurations have the low- 
est energy. This is in strong contrast with what hap- 
pens in the hole doped case where inhomogeneities are 
already found very close to half-filling22,. The result is 
compatible with the recent theoretical proposal of phase 
separation at intermediate fiUings^^:. Phase separation 
has also been claimed for the electron doped t — t' — J 
modcl^. Unrestricted Hartree-Fock can work better to 
explore the electron doped cuprates than the hole doped 
cuprates because when doping the t — t' Hubbard model 
with electrons we fill the wider upper antiferromagnetic 
band where quantum fluctuations are expected to be less 
important. 
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ARPES experiments in electron doped^'' cuprates show 
a Fermi surface made of electron pockets for dopings of 
X = 0.04, X — 0.10 and x — 0.15 as expected from 
long range antifcrromagnctic orderiS^. They also see the 
emergence of a hole-like Fermi surface around the nodal 
direction at a; = 0.15. This nodal signal seems to indi- 
cate the instability of the antiferromagnetismiSiSl. As 
we work in the strong coupling limit we do not obtain 
the hole-like Fermi surface at a; = 0.15. It is known that 
to get it in mean-field calculations unrealistic values of 
U are neededS^. This case has been explored with dif- 
ferent techniques where quantum fluctuations are prop- 
erly taken into account (seei^ii^ii^) but these techniques 
are not well suited for studying inhomogeneities. On the 
other hand ARPES also indicates the presence of spectral 
weight in the insulator gap. Theoretically, these midgap 
states have been obtained by adding spin fluctuations to 
the mean field solution^^. Inhomogeneous configurations 
also form midgap states which would be an alternative 
explanation of these states. 



A related issue concerns recent puzzling experimen- 
tal results of electron-doped cuprates. Although the 
Fermi surface found in ARPES for most of the electron- 
doped cuprates lies in the antinodal direction, recent 
Raman experiments have found that coherent quasi- 
particles mainly reside in the nodal direction of the Bril- 
louin zona^^ in all the samples that they have analyzed 
{x = 0.135, X = 0.147 and overdoped samples). In partic- 
ular they see that the B2g signal which involves averages 
mostly over the nodal region has a stronger quasiparti- 
cle component than the Big signal which deals with the 
antinodal. Both signals show an incoherent background. 
In our calculations we find cases of inhomogeneous config- 
urations that present a quasiparticle-like signal in the i?2g 
Raman response although the Fermi level does not cross 
the nodal region. We argue that this is an example of a 
case in which two-particle properties cannot be induced 
from one-particle properties. We also find that inhomo- 
geneities provide an incoherent background to the Raman 
response absent in the case of homogeneous metallic an- 
tiferromagnetism. Our results are qualitative. We do 
not pretend to do a close comparison with experiments 
for what we would have to take into account other very 
important effects as quantum fluctuations, scattering by 
impurities etc. We intend to find out the prints of inho- 
mogeneities in the Raman signal. 



II. THE MODEL AND THE METHOD 

The t — t' Hubbard model is defined in the two dimen- 
sional squared lattice by the hamiltonian 

H = -t c+cjs -t' ^ c+cjs + U^n,^nii, 

<i,j>s <<i.,j>>s i 

(1) 

which gives rise to the dispersion relation 

e(k) ~ — 2t [cos{kxa) + cos{kya)] — At' cos{kxa) cos(fcj,a). 

t > and t' = —Q.3t are generic values for modeling the 
physics of the cuprates. We will consider electron doping 
{n > 1). We work in lattices of different sizes with pe- 
riodic boundary conditions depending on the problem at 
hand. We obtain the U — n phase-diagram working in a 
L = 14 X 14 lattice. We use a 16 x 16 lattice to discuss in- 
phase versus off-phase configurations and a L = 24 x 24 
lattice to compute the spectral function and the Raman 
response. A comparison of the results on the different 
lattice sizes for homogeneous and inhomogeneous config- 
urations is done in the next section. The unrestricted 
Hartree-Fock approximation minimizes the expectation 
value of the hamiltonian in the space of Slater de- 
terminants. These are ground states of a single particle 
many-body system in a potential defined by the electron 
occupancy of each site. This potential is determined self- 
consistently 

H — — y ^ tjjcl^Cjs — y ] —friicl^ass'Cis' 

i,j,s i,s,s' 

+C.C., (2) 

i 

where tij denotes next, t, and next-nearest, t', neighbors, 
and the self-consistency conditions are 

TOi ^ < cl^ass'Cis' >, qi =< n^i + - 1 >, 

where a are the Pauli matrices. 

We will call E\ where A = 1 . . . 2_L the eigenvalues 
of the hamiltonian and u||0), . . . ,U2^|0) the basis that 
diagonalizes the hamiltonian. These fermionic operators 
are obtained from the original ones through: 

Cjs =^Ujs,XUx. (3) 
A 



The organization of this paper is as follows. In section 
II we review the t — t' Hubbard model and the unre- 
stricted Hartree-Fock method. In Sect. Ill we discuss 
the U — n phase diagram. Sect. IV gives the results 
obtained on one-particle properties and Section V deals 
with the electronic Raman response. In section VI we 
present our conclusions and open problems. 



III. U -n PHASE DIAGRAM 

We have built a detailed U — n phase diagram for 
t' — — 0.3t, which it is shown in Fig. ^ by comparing 
the energies of paramagnetic, antiferromagnetic and in- 
homogeneous configurations in a 14 x 14 lattice size. An- 
tiferromagnetism prevails over most of the phase diagram 
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FIG. 1: U — n phase diagram. The symbols represent the 
values studied. AF stands for antiferromagnetism, VS for in- 
phase vertical stripes, the inhomogeneous configuration rep- 
resented in the Fig. |3|(up), and PM for paramagnetism. At 
strong coupling and intermediate doping inhomogeneous con- 
figurations (VS) are the dominant ones in the phase diagram. 



and inhomogeneities (VS in Fig. ^ ) can be found at in- 
termediate doping and for high values of U. For the in- 
homogeneous solutions we have compared diagonal (DS) 
and vertical (VS) in-phase stripes, in-phase antiferromag- 
netic domains (DOM) and polarons (POL). Fig. [21 shows 
the comparison of the energies per particle for U — 8t 
and t' = —0.3t. All these inhomogeneous solutions con- 
verge for n = 1.15 where no antiferromagnetism can be 
found. In fact, antiferromagnetism does not converge in 
all the VS region in the phase diagram of Fig. ^ This 
is a clear indication that there is a critical density where 
inhomogeneities are preferred over homogeneous antifer- 
romagnetism in unrestricted Hartree-Fock. The inhomo- 
geneous solution lowest in energy consists of in-phase ver- 
tical stripes (VS) represented in Fig. |21(up) in a 24 x 24 
lattice. In-phase antiferromagnetic domains (DOM) are 
very close in energy and converge very easily. They are 
represented in Fig. |3| (down). Magnetic polarons are only 
stable in the point U = 8t and n = 1.15. Away from it 
they are such that even if we start to run the system in a 
polaronic configuration it ends up in the metallic antifer- 
romagnetism, i.e polarons dilute in the antiferromagnetic 
background. This is in contrast with the case for t' — 
where polarons are the preferred configuration^^ at low 
and intermediate dopings for high values of U. The in- 
phase diagonal stripe (DS) only converges at U — 8t, 
t' ~ -~0.3t and n — 1.15 but it is the highest in energy, 
however DS is the lowest in energy in hole doped Hub- 
bard model within UHF^"'^^. We will further comment 
on these results when we discuss the density of states. 

To check that the results obtained are robust over the 
size of the lattice and as in the following section we will 
use results for a 24x24 lattice we have made a comparison 
of the energy per particle for different lattice sizes. The 
results are shown in Table ^ The compared energy per 
particle are for homogeneous antiferromagnetism (AF) 
and inhomogeneous in-phase stripes (VS) configurations 




1 1.05 1.1 1.15 1.2 1.25 



n 

FIG. 2: (Color online) Comparison of the energy per particle 
for different configurations in a 14 x 14 lattice for different den- 
sities n. AF stands for antiferromagnetism, VS for in-phase 
vertical stripes (Fig. El (up)), DOM for in-phase antiferromag- 
netic domains (Fig. |2l (down)), POL for polarons and DS for 
in-phase diagonal stripes. The inset shows the U — St and 
n = 1.15 region enlarged to appreciate better the difference 
in the energies. 



Conf. E/L{12 x 12) E/L{2A x 24) 
AF -4.46587799 -4.46587795 
VS -4.04248799 -4.04257066 

TABLE I: Energy per site in units of t of different configura- 
tions for different lattice sizes 12 x 12 and 24 x 24 for U = 8t 
and t' — —0.3t. In the homogeneous solution (AF) at n = 1 
the difference is negligible and for the inhomogeneous solution 
(VS) at n = 1.16 it is appreciable in the fourth digit. 



for two lattice sizes 12 x 12 and 24 x 24. The lattice size 
has been chosen for convenience to fit the inhomogeneous 
configuration. In the homogeneous case the variation in 
energy is almost negligible and in the VS case it is ap- 
preciable in the fourth decimal digit. Therefore we can 
trust that the phase diagram will remain for bigger lattice 
sizes. 

Next we compare in-phase with off-phase stripes. We 
have chosen a different lattice size (16 x 16) to fit the 
off-phase configuration. For the points shown in Fig. ^ 
off-phase vertical stripes only converge at n = 1.15 and 
U = 8t as an excited state, though the energy per particle 
is very similar to the in-phase stripe {Eys/L = — 4.090i, 
EsTRiPE I L = — 4.083i). Therefore, within Unrestricted 
Hartree Fock, in-phase stripes are preferred. Adding fiuc- 
tuations to the unrestricted Hartree-Fock solution could 
change this view. Calculations using an exact diagonal- 
ization method within the dynamical mean field theorji^ 
have shown that solitonic formation contributes to the 
stability of the off-phase stripes. However we have no- 
ticed a clear tendency for in-phase solutions that con- 
verge easily and in a wider region of the phase diagram. 
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FIG. 3: Configuration for U = 8i, t' = -0.3t and n = 1.15 in 
a 24 X 24 lattice for in-phase vertical stripes -VS- (up) and in- 
phase antiferromagnetic domains -DOM- (down). The circles 
and arrows are scaled to the local charge and spin. 



IV. ONE-PARTICLE PROPERTIES 

The spectral function ^(k, cj) = — i3Gret(k, w) is 
evaluated using the time-dependent Green's function, 
which can be calculated numerically. When the spec- 
tral function is expressed in terms of the new manybody 
state by Eq. Othe following expression is obtained" 

m.n,s.A 

(4) 

We represent A(k, uj) in Fig. 0] with the Fermi level 
at w = 0, for AF (Fig. Hup) at n = 1.05 and for DOM 
(Fig.^down) at n = 1.15. The spin-charge configuration 
of DOM is represented in Fig. |31 (down) in a 24 x 24 lattice 
for parameters U = 8t and t' = —0.3t. The VS repre- 
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FIG. 4: (Color online) A(k, ui) for different configurations an- 
tiferromagnetism at n = 1.05 (up) and DOM at n = 1.15 
(down). Parameters are U = 8t, t' = -0.3t and T = 0.05t. 

sented in Fig. (up) has a similar spectral density as the 
DOM. We represent the DOM spectral function instead 
the VS because this result will give us a coherent quasi- 
particle component in the i?2g Raman response though 
the Fermi level does not crosses the band in the nodal 
region. In Fig. ^a) we observe two well known results, 
(l) t' makes the lower band narrower and the upper band 
wider. In consequence, it is expected that localization ef- 
fects become more important when we dope with holes 
than when we dope with electrons. (2) concerning the 
k dependence for homogeneous antiferromagnetism, at 
low doping the Fermi level goes to the bottom of the up- 
per band at (tt, 0) as suggested in a rigid band structure. 
This result is compatible with ARPES2&. For inhomo- 
geneous configurations. Fig. 0fb) , we do not longer have 
the rigid band structure view and the spectral weight 
is reorganized. Midgap states appear and the spectral 
weight around (tt, tt) in the upper band disappear. The 
Fermi level is again at the bottom of the upper band at 
(tt, 0). The midgap states will be more appreciated in the 
density of states. 

Fig-Elshows the density of states for the AF and DOM 
configurations described previously and also for the VS 
represented in Fig.|21(up) with a width of 77 = O.lt given 
to the delta functions. The Fermi level is located at zero 
energy. We again see for the antiferromagnetic solution 
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FIG. 5: Density of states for homogeneous antiferromag- 
netism at n = 1.05 (dotted line), in-phase antiferromagnetic 
domains DOM (solid line) and in-phase vertical stripes VS 
(dashed line) at n = 1.15. Both inhomogeneous configura- 
tions present similar density of states with in-gap states and 
a pseudogap-like feature at the Fermi level. Parameters are 
U = St, t' = -0.3t, r = 0.05t and r; = O.lt. 



(dotted line) the result expected from a rigid band struc- 
ture, the gap is lower if we compare with the result at 
half filling (not shown) and the Fermi level is shifted to 
the upper band. For DOM (solid line) and VS (dashed 
line) the density of states is very similar. We appreci- 
ate the formation of midgap states with almost no gap 
closing if we compare with the former case at lower dop- 
ing 71 = 1.05. It is interesting to notice in Fig. [S] that 
for the inhomogeneous solution there is a depression at 
the Fermi level (it can also be seen in Fig. Q^down)). 
We do not claim that this corresponds to the pseudogap 
whose physics needs better approaches, but it is inter- 
esting to realize that inhomogeneous solutions can give a 
feature that it is usually considered as a hallmark of the 
pseudogap>i& 

Next we discuss further the phase diagram of Fig. 
with the help of the spectral function (Fig. ^ and the 
density of states (Fig. [SJ . A simple argument to under- 
stand why the antiferromagnetism is enhanced and dom- 
inates the phase diagram of the electron doped case can 
be given by changing to a Neels state basis. In that basis 
the t' term appears in the diagonal of the Hamiltonian 
since it represents a hopping within the same sublattice. 
As such it does not change the antiferromagnetic back- 
ground. The first electron carrier will go to (tt, 0) (see 
Fig. EI up)) with an energy given by i?(7r, 0) — 4t' + A 
where A = mil and m is the magnetization. As t' < 
it lowers the energy of the antiferromagnetic solution. A 
similar argument is given in Refi^ for the electron doped 
t — J model. This simple argument is reinforced by the 
study of the stability of the metallic antiferromagnetic 
solution^i. In the density of states represented in Fig. [S] 
can be seen that our antiferromagnetism is also metallic 
and the whole picture is very similar to what happens in 
a rigid band model. 

Next we discuss the VS at intermediate doping and 



high U. As mentioned when doping with holes we are 
doping a narrower band and a lower mobility is expected 
than when doping with electrons. It is not surprising 
then that inhomogeneities appear at low doping in hole 
doped and not in the electron doped. However there is a 
critical density in the electron doped case where inhomo- 
geneities are the preferred configuration as a compromise 
to minimize the strong Coulomb interaction. In-phase 
vertical stripes allow more mobility than diagonal stripes 
and that is why within UHF diagonal are preferred for 
hole doped^SiSi and vertical for electron doped. 



V. NON-RESONANT RAMAN RESPONSE 

The differential Raman cross section is^ 



(5) 



where rg = /mc^ is the classical radius of the electron 
and the generalized dynamic structure factor is 



5^^(q,w,/3) = ^ ■ 
hf 



m~p^(Si)\i)\H{Ef-E, + h^). 



Here, 



(6) 

— is the probability of the initial state in ther- 
modynamic equilibrium. The state \i) is the initial many- 
body state and |/) is the final many-body state. p(q) is 
an effective charge density given by 



(7) 



where s is the spin index and 7(k; w/, ws) is the Raman 
scattering amplitude. We are interested in the q = 
case since the momentum transfer by the photon is neg- 
ligible. The generalized structure factor given by Eq. |^ 
can be interpreted as the power spectrum of the charge 
density fluctuations at frequency uj. It is related to the 
imaginary part of the Raman response function through 
the fluctuation-dissipation theorem 



where 



^^.^(q,^) 



X77(q,cj) = 



-[1 + n(tj)]Imx77(q, w), 



l/T 



(ire 



■(r.p(r)p(O)). 



(8) 



(9) 



We will compute the imaginary part of the Raman re- 
sponse by replacing the Hartree-Fock eigenvalues and 
eigenvectors in the generalized structure factor given by 
Eq. 1^1 instead of using the Green's function machinery in 
Eq. |n|as it is usually done. 

We will follow Ref4£ where symmetry is used to clas- 
sify the scattering amplitude since the different polar- 
izations transform as elements of the point group of the 
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crystal 



7(k;w/,tJs) 



L 



(10) 



where $L(k) are the Brillouin zone harmonics (BZH). 
The use of symmetry arguments instead of the inverse 
effective mass approximation as it is usually dona^S will 
allow us to discuss non-resonant Raman response at high 
energies^. In this work we are interested in the B^g and 
Big symmetries of the square lattice. The Big symmetry 
corresponds to incident and scattered light polarizations 
aligned along x -\- y and x — y, and thus 



$Big(k) = cos(kx) - cos(ky 



(11) 



where the dots stand for higher order BZH. The Big sym- 
metry corresponds to incident and scattered light polar- 
izations aligned along x and y 



^B2g (k) = sin(kx) sin(ky) 



(12) 



The Big signal weights mainly the antinodal quasiparti- 
cles and the Big signal the nodal quasiparticles. 

Replacing Eq. O in Eq. [71 and Eq. El we obtain for the 
Big and the Big Raman expressions 



-l3Ex 



, (l + e-/5^v)(l-He-/'^^) 



5{io + {E^-E^,))\ J2 iUl+i,y)s.,U^.,y)s,x' + h.c) 

S = ix.y)s 

-iUl^y+i)s,xUi..y)s,y + h-c)\\ (13) 



\ (l + e-'3i^A')(l + , 



) 



h.c) 



Fig. ini shows the non-resonant Raman response. In ex- 
periments for cuprates, the range in frequencies is just 
up to around It and at high frequencies the resonant Ra- 
man response is unavoidable and it could couple with 
the non resonant one as has been recently shown with 
Dynamical mean field theory*^, but we show the com- 
plete range of frequencies to illustrate the physics that 
comes up from our calculations. We can phenomeno- 
logically decompose the Raman response as following: 
X"(w) = Xio^^ + Xincoh + Xgap- The low frequency 
response is usually fitted as a quasiparticle response 
XQEpi'^) ~ '^ ui'^+r'^ ' '^^^ r = since we did not 
put any scattering rate as input, and therefore Xqep — ^■ 
But for some cases we find xiowuj 7^ as we will see in 
the following, x'gap formed by the pronounced peak 




FIG. 6: (a) Evolution of the Big and B2g Raman signal with 
doping for antiferromagnetism for U = 6t, x = 0.05 (solid 
line) and a; = 0.1 (dashed line), (b) Big and -B29 Raman 
signal for VS (Fig. E] (up)) at [/ = 8t and n = 1.15. (c) 
Big and -623 Raman signal for DOM (Fig. (down)) at [/ = 
8t and n = 1.15. In the inset is shown the evolution with 
temperature of the signal for low frequencies. Parameters 
are t' = -0.3t, T = 0.05t and r] = O.lt. 



at high energies in both figures and is due to the an- 
tiferromagnetic gap. Xincoh formed by the rest. In 
experiments this part is almost a constant continuum 
background and gave rise to the Marginal Fermi liquid 
phenomenology^^. Next we discuss the results. Fig.jSJa) 
represents the evolution with doping of the Big and Big 
Raman signals for homogeneous antiferromagnetism with 
parameters U — 6t, t' = —Q.3t and dopings n = 1.05 
(solid line) and n = 1.1 (dashed line). Again a width 
of J] = O.lt is given to the delta functions. As expected 
for homogeneous antiferromagnetism the signal for the 
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channel is negligible and cannot be appreciated in 
Fig. EJa) and in the Big channel, Xgap is the only ap- 
preciable component. We can observe a shift in Xgap for 
n = 1.1 compared with n = 1.05 reflecting the closing 
of the gap with increasing doping. That x'l'owu: ^^"^ 
x'incoh ^'^^ negligible might be due to the fact that the 
charge is completely homogeneous for this antiferromag- 
netism and from Eq. ^| and Eq. E| we see that for ho- 
mogeneous phases is very likely that the Raman matrix 
elements are zero. We have checked this result with all 
the homogeneous configurations that we have found with 
the same nuU resuh for x'l'owc, and x'Lcoh- Fig-EIb) and 
(c) represents the Big and B2g Raman signals for the 
inhomogeneous solution consisting of VS and DOM at 
U = 8t, t' = —0.3t and n = 1.15. In both figures a gap is 
observed at high frequencies and x'incoh longer zero 
neither in Big nor in B2g channels. The gap appears at 
higher energy in the Big channel consistent with Fig. ^ 
(down). The incoherent background comes from the in- 
gap states in the density of states, see Fig. El For VS 
there is not low-frequency component in any channel but 
for DOM, the i?2g channel has a x'lowuj component with 
a quasiparticle-like shape and it is not clear in the Big 
channel. We have calculated the evolution with temper- 
ature of the Raman signal shown in the inset of Fig.EJc). 
The decrease with temperature of the Xqep ^o'' ^'^g 
channel would have been interpreted as a metallic behav- 
ior in experiments. We arrive at two interesting results 
for the inhomogeneous solution, (1) there is a remarkable 
incoherent component, and (2) the B2g channel has a 
quasiparticle-like component for the DOM configuration. 
We can also find this behavior for higher doping in the 
VS configuration. This means that the low energy part 
of the Raman response can arise not only from the real 
quasiparticle component but also from inhomogeneities. 
In the last case this quasiparticle-like response has no re- 
lation with real quasiparticles since the Fermi level does 
not cross the nodal region. It might arise because in- 
homogeneities create an effective Raman scattering rate 
that can contribute to the low frequency region of the 
Raman response. 

VI. CONCLUSIONS 

We find in-phase inhomogeneities as the lowest en- 
ergy configurations for strong-coupling limit at interme- 



diate dopings. These results are interesting because the 
doping region where inhomogeneities are found signals 
another quantitative asymmetry for hole and electron 
doped cuprates. Inhomogeneities are present at very low 
doping for hole doped cuprates and only at intermedi- 
ate to high doping for electron doped. This region of 
intermediate doping is precisely the region where the an- 
tiferromagnetic phase is suppressed {n « 1.5). The result 
agrees with a recent theoretical proposal of phase sepa- 
ration for this range of parameters'^. It is also interest- 
ing to notice that we need to be in the strong coupling 
limit to find the inhomogeneities, since they are formed to 
minimize the strong on-site interaction. Experimentally 
electron doped cuprates seem to be in the strong cou- 
pling limit as much as hole doped cuprates^-"'^ but there 
is a theoretical debate on whether or not strong coupling 
is necessary to understand electron doped^^. In-phase 
inhomogeneities are compatible with experiments of in- 
clastic neutron scattering and thermal conductivity^^. 
However to fully address this issue quantum fluctuations 
should be taken into account since they can favor off- 
phase stripes^. 

The inhomogeneities give rise to midgap states in the 
density of states as found by ARPES^^ though at strong 
coupling it is not possible to obtain a quasiparticle signal 
in the nodal region, for what better techniques taking 
into account quantum fluctuations should be used. For 
the non-resonant Raman response the in-phase inhomo- 
geneities give an incoherent background at low and inter- 
mediate frequencies. For the DOM configuration (Fig. |3| 
(down)) the i?2g signal have stronger quasiparticle-like 
component than the Big signal though the Fermi level 
does not cross the nodal region. This is an example where 
two particle properties can not be inferred from one par- 
ticle properties. It would be interesting to check if the 
quasiparticle component of the i?2g Raman response is 
different from zero in the region where ARPES does not 
show quasiparticle weight in the nodal region. 
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